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In this paper, we thoroughly investigate the LHC phenomenology of the type II seesaw mech¬ 
anism for neutrino masses in the nondegenerate case where the triplet scalars of various charge 
, H^, , A°) have different masses. Compared with the degenerate case, the cascade 

decays of scalars lead to many new, interesting signal channels. In the positive scenario where 
Mh±± < Mfj± < AI[jo/ao, the four-lepton signal is still the most promising discovery channel 
for the doubly-charged scalars The five-lepton signal is crucial to probe the mass spectrum 

of the scalars, for which, for example, a 5a reach at 14 TeV LHC for Mfj± = 430 GeV with 
Mh±± = 400 GeV requires an integrated luminosity of 76 fb^^. And the six-lepton signal can be 
used to probe the neutral scalars /A^, which are usually hard to detect in the degenerate case. In 
the negative scenario where Mh±± > Mfj± > M^o /^o, the detection of is more challenging, 
when the cascade decay H^W^* is dominant. The most important channel is the associ¬ 
ated /A^ production in the final state C^^Tbbhh, which requires a luminosity of 109 fb~^ for 

a 5cr discovery, while the final state C^^TbbT'^T~ is less promising. Moreover, the associated H'^A^ 
production can give same signals as the standard model Higgs pair production. With a much larger 
cross section, the H^A^ production in the final state bbT~^T~ could reach 3a significance at 14 TeV 
LHC with a luminosity of 300 fb“^. In summary, with an integrated luminosity ~ 0(500 fb~^), the 
triplet scalars can be fully reconstructed at 14 TeV LHC in the negative scenario. 
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I. INTRODUCTION 


When the standard model (SM) is considered as a low energy effective field theory, the tiny neutrino 
masses can be incorporated by higher dimensional operators that are suppressed by new physics above 
the electroweak scale. Such an operator first appears at dimension five, which is fhe unique Weinberg 
operator lli, O 5 = (F^ed>) where Fl and are respecfively fhe leff-handed lepfonic and Higgs 

doublefs in SM and e is an anfisymmefric mafrix. If fhe operator is suppressed for one reason or anofher, 
fhe neufrino masses would be induced by even higher dimensional operafors, which are also known fo be 
unique af each dimension, 05 + 2 n = wifh n aposifive integer {21. 

If is inferesfing fhaf fhere are exacfly fhree possible ways af free level to realize fhe Weinberg opera- 
for in an underlying fheory Q- They correspond fo fhe fype I JH, II @, and III f6] seesaw mechanism 
fhaf infroduces respecfively fermionic singlefs, a scalar friplef of hypercharge 2, and fermionic friplefs of 
hypercharge zero. Varianfs of fhe fhree seesaws are also suggesfed based on differenf fheorefical consider¬ 
ations Q. For insfance, a combinafion of fype I and III seesaws has a simple realizafion in fhe confexf of 
grand unified fheories fT], while leff-righf symmefric models |'81 embrace fype I and II nafurally. The liny 
mass of neulrinos can be more readily induced if fhe effecfive operafors are radiafively induced f9] or firsl 
appear af free level af a higher dimension |[T0l . 

To discern fhe underlying physics fhaf is responsible for neufrino masses, if is vifal to produce direcfly 
the relevant heavy particles at colliders and study their properties |[TTI - [2^ . In this paper, we will study the 
LHC phenomenology of the type II seesaw mechanism. We recall in this circumstance that the previous 
literature mainly concentrates on the simplified version of if, assuming fhaf fhe friplef scalars are degenerafe, 
and on fhe search of fhe doubly charged scalars ifT^ IT4l [TTl |25| - IT7]| . Assuming degeneracy, 
have only Iwo decay modes: fhe leplon number violafing (LNV) like-sign dileplon decays 
and fhe like-sign diboson decays —)• W^W^. Given fhe range of neufrino masses, Iheir relalive 

imporlance is confrolled by fhe vacuum expecfalion value (vev) of fhe friplef, V/\. The dilepfon decays 
dominafe when ua < 10“^ GeV, resulfing in fhe characferisfic four-lepfon signal from pair production 
pp —)■ —)■ ifT^ [TTl l26l and also fhe promising fhree-lepfon signal from associated 

production pp —)• —)■ llT4l[T7ll26l . For va > 10“'^ GeV instead, fhe diboson decays are 

dominanf, and fhe imporfanf signal channel is pp —)• H~^^H —> W^W~^W~W~ —)■ ifT^ 

as well as associated production pp —)■ —)• IHl- 


The (dominanfly) friplef scalars are generally nondegenerafe, and fo good accuracy fheir mass spliffings 
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are determined by a coupling constant A 5 in the potential, see Eq. Q. We discriminate two scenarios: 
positive scenario (A 5 > 0) : Mh±± < Mii± < 

negative scenario (A 5 < 0) : Mfj±± > Mu± > (1) 


The collider signatures for the nondegenerate case have received relatively less attention. For the positive 
scenario, Ref. computed the enhancement of pair production due to cascade decays of 

heavier scalars and mentioned further contributions from cascade chain decays of H^/A^ production. 

The cascade chain decays can produce a like-sign four-lepton signal whose production rate was estimated 
in Refs. 12^ 1^ in the narrow width approximation. The authors in Ref. |[^ investigated the five- and six- 
lepton signals from cascade decays of /A^, and H^A^. For the negative scenario. Ref. |[32l 

noticed the importance of the cascade decay H^W^* and compared its branching ratio with 

dilepton decays —)■ The work ll^ focused on the signal channels pp ——>• £~^^Tbbbb 

for the search and pp —)• £~^i^^Tjjbbbb for the search; see 13^ for more details. 

Nevertheless, to the best of our knowledge, a systematic and comprehensive study of FHC signatures is still 
lacking. The purpose of this paper is to fill fhe gap. 

The resf of fhe paper is organized as follows. In Sec. we sef up our nofafions in fhe fype-II seesaw 
model, and review fhe currenf experimenfal consfrainfs on fhe model paramefers. The decay properfies 


of fhe friplef scalars in fhe nondegenerafe case are invesfigafed in Sec. Ill This is fhen followed by fhe 
core Sec. IV and Sec. |Vj where we sfudy sysfemafically fhe FHC signafures for fhe positive and negafive 
scenarios in fhe nondegenerafe case. We esfimafe fhe SM backgrounds and develop fhe sfrafegies fo separafe 
fhe signals from backgrounds in each signal channel. Finally, in Sec. |V^ we presenf our conclusions wifh 
critical discussions. 


II. CONSTRAINTS ON PARAMETERS IN TYPE II SEESAW 


We will review in fhis secfion fhe currenf experimenfal consfrainfs on fhe fype II seesaw model, so fhaf 
our lafer phenomenological analysis af FHC can be more realisfic. To sef fhe sfage, we firsf give a concise 
infroducfion fo fhe model. The seesaw operafes wifh fhe help of a scalar friplef A of hypercharge 2 in 
addition fo fhe SM scalar doublel of hypercharge 1: 


/ (/)+ \ ^ / 5 +/V 2 ,5++ \ 

V ’ V 5 ° -S+/V2) 


( 2 ) 
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where the superscripts denote the electric charge. The most general potential is given by 

F(4>, A) = + M2Tr(AtA) + + Aa (Tr(AtA))^ + A3Tr(AtA)2 

+A4(4>^4>)Tr(AtA) + A54>'^AAtT>+ (^^4>'^ir2A+4> + h.c.) , (3) 


where spontaneous symmetry breaking is triggered by assuming < 0 that results in a vev (v) for the 
scalar doublet, while the parameter is assumed to be positive to set the mass scale for the heavy scalars. 
The /i term is important, and deserves a few comments. It induces a vev (ua) for the triplet out of that 
for the doublet. This in turn causes mixing between the doublet and triplet scalars of equal charge. Since 
the /r parameter can be taken real without loss of generality, CP is preserved by the potential and mixing. 
Together with the Yukawa couplings between the leptons and the scalar triplet, it violates the lepton number 
by two units. The fj, parameter and thus va are considered to be naturally small in this sense. The triplet 
vev causes the deviation of the p parameter from unity at tree level, i.e., p k, 1 — 2v\jv^ for va ^ v. 
The precise experimental measurement on p then translates to a limit on ua, which we take safely to be 
UA < 1 GeV in our numerical analysis. 

Separating out the vev’s. 


+ (/> + ix), 5° = + S + iO, 


the scalars mix as follows (see for instance. 


(4) 


4>^\ / 


X \ /GO 

' =Ria) 

u V 


h 

= RiOo) , ^ 

6 1 \ 


(5) 


Here R{u}) is the standard rotation matrix in the plane, and the mixing angles are given by 

^ a V^VA ^ 2 va ^ VA 2u 2(A4 + As) - 4M| 

tan6»+ =-, tana =-, tan26'o = —A 2 \' " ' 2 A i x V' 

V V V 2v^Xi - - v^[X2 + X 3 ) 


( 6 ) 


The physical scalars thus include the doubly-charged H^^{= <5^^), singly charged and the CP-even 
(-odd) neutral h, (^O), while are the would-be Goldstone bosons. An auxiliary parameter is 
introduced for convenience 


Af 2 _ 

^ x/ 2 ua^ 


(V) 


so that the heavy scalars have the masses approximately 


M^±± ^Ml- ^A5^;^ M^± « Mi - ^A5r;i M|o ^ Mi, 


( 8 ) 
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while the light scalar has the mass M| « 2Xiv^. The heavy scalars are equidistant in masses squared to 
good approximation: 

M^±± - Ml± « M^± - M|o (9) 

with A^, and therefore we distinguish between the two scenarios of masses as shown in Eq. 

Q according to the sign of A 5 while the degenerate spectrum corresponds to the special case of A 5 = 0 . 
There are various theoretical considerations such as perturbativity, vacuum stability, and unitarity that are 
employed to constrain the parameters in the potential. A detailed study of the potential V (<h, A) shows both 
scenarios are allowed (351. The parameters chosen for our collider simulation will be located in the allowed 
region. 

Since the triplet scalars contribute through radiative corrections to the precisely measured electroweak 
quantities, their parameters are constrained by the electroweak precision data. A careful analysis of one- 
loop radiative corrections has been made in Ref. (3^ . and the result turns out to depend significantly on 
the renormalization scheme employed. In the so-called scheme I, one employs the effective mixing angle 
defined by the Zee vertex as the fourth input parameter in the gauge sector. This is the scheme also 
used in Refs. (3411491 . It was found that the measured value of mw requires a large mass splitting AM in 
the positive scenario for ua ^ 1 GeV, Mjj++ ~ 150 — 300 GeV and specific values of ofher paramefers, 
while fhere is severe tension in fhe negafive scenario for fhe same ranges of parameters. Here AM is 
defined as fhe (positive) mass splitting befween and fhe lighfer of and In fhe scheme II, one 
adopfs instead fhe mixing angle a befween and as fhe fourfh inpuf paramefer. If was found (36l 
fhaf for UA ^ 1 GeV fhe posifive and negafive scenarios can be accommodafed wifh AM < 50, 30 GeV 
respectively. This resulf is consisfenf wifh fhe bound from fhe elecfroweak S, T, U paramefers (37l . which 
requires AM < 40 GeV, independenfly of fhe doubly charged scalar mass. As we are interested here in 
fhe case of small ua, we will assume Ibis latter sef of consfrainfs on AM in our physics analysis. 

The Yukawa coupling befween fhe scalar friplef and lepfon doublefs is responsible for neufrino masses: 

^Yuk = -YijFu + h.c., ( 10 ) 

where fhe superscripf C denofes charge conjugafion and i, j are generafion indices. The mafrix Y is 
generally complex and symmefric, and gives fhe Majorana neufrino mass mafrix 

M^ = V2Yva. (11) 

In fhe basis where fhe mass mafrix of charged lepfons is diagonal, fhe neufrino mass mafrix is digonalized by 
the unitary PMNS matrix, Mjy = fpMNs"^i^^MNS’ being diagonal, real, and semi-definite positive. 
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The matrix Y is thus governed by the neutrino spectrum, mixing pattern, and the triplet vev va, and the 
relationship is constrained through Yukawa contributions to the leptonic decays of the charged scalars and 
the lepton flavor violating (LFV) transitions of the charged leptons. 

The mixing matrix Vrmns is generically parameterized in terms of three mixing angles 0ij, a Dirac CP 
phase 6 , and two Majorana CP phases. Assuming vanishing Majorana phases, we use the following best fit 
values lIMll of parameters for the normal (inverted in parentheses, if different) hierarchy of masses: 

Amii = 7.62 X 10"'^ eV^ \Aml^\ = 2.55 (2.43) x 10"^ eV^; 

sin^ 012 = 0.320, sin^ 023 = 0.613 (0.600), sin^ 0i3 = 0.0246 (0.0250); 

6 = O.SOvr (-0.037r). (12) 


The absolute neutrino masses remain unknown. Cosmological considerations set an upper bound on the 
sum of masses 1391 |40l, Yli < 0.23 eV. The neutrinoless double-/? decay is sensitive to the effective 
Majorana mass IHTII . (m)ee = | most stringent bound coming from the EXO 


Colla. i42l . (m)ee < 0.14 — 0.38 eV. The direct neutrino mass search is based on kinematics and sensitive 
to the average electron neutrino mass, = Yli I (^Ns)eiP"ii ■ The current record is kept by the Troitsk 
Colla., rriu^ < 2.0 eV |43| . and the upcoming KATRIN experiment has the potential to reach a level of 
0.2 eV i44|. Considering all of these constraints, we assume the lightest neutrino to be massless in either 
normal (NH) or inverted hierarchy (IH). 

The Yukawa couplings can mediate low energy LFV processes. The purely leptonic decays 4 —)• ijik^i 
proceed at tree level by the exchange of the doubly-charged scalars while the radiative transitions 

(-jl gain contributions from virtual and loops; see, for instance. Ref. Il45l for a compre¬ 
hensive analysis. The parameter space of My to be probed with such processes in ongoing and planned 
experiments has been discussed in Ref. Il46ll . Using the experimental limits on LFV processes and the muon 
anomalous magnetic dipole moment and assuming degenerate heavy scalars for simplicity. Ref. WT\ gives 
the combined bound, vaMu±± > 150 eV GeV. Since both and contribute to, e.g., /r —)■ ey, 

|2 r 

I I 

(13) 


\ «QEd|(M, 


F^A 


M‘^±± 


+ 


8M|±j 


the bound will be different in the nondegenerate case. Using the current most stringent bound by the MEG 
Colla., BR(^ —)• ey) < 5.7 x 10“^^ (90% C.L.) 1481 . we draw on the left panel of Fig.j^the lower bound on 
VA as a function of Mh±± for a given mass splitting of the charged scalars 6 M = Mh± — Mh±± for both 
NH and IH. {6M = AM in the positive scenario and \6M\ ss AM in the negative scenario for a relatively 
small mass splitting.) One sees that the bound for IH is generally more stringent than for NH, and that the 







7 


positive (negative) scenario is more loosely (stringently) constrained than the degenerate case: 


6 M = 40 GeV 
6 M = 0 GeV 
5M = -40 GeV 


vaMh±± > 180(NH),235(IH) eV GeV; 
vaMh±± > 185(NH),240(IH) eV GeV; 
vaMh±± > 200(NH),260(IH) eV GeV. 


(14) 


On the right panel of Fig. we show the relative change by the ratio Rsm = va{SM) /va{0)- The deviation 
between the positive scenario and degenerate case is small and usually cannot excess 3%, while it can reach 
about 10% for the negative scenario. For heavier than 200 GeV, all deviations can be safely neglected 
for both scenarios. Finally, we should mention that the above limits also depend on the lightest neutrino 
mass mi (NH) or m 3 (IH) and Majorana phases that we have assumed to be zero. 
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FIG. 1. Lower bound on va for nondegenerate case (left panel) and ratio Rsm (right) as a function of M^±±. 


We finally turn to the constraints at colliders. The decay /i —)■ 77 played a central role in the discovery 
of h. While the tree-level decays of h keep essentially intact in the type II seesaw due to small mixing 
governed by va ^ v, the loop-induced decays h —)• 77 , Z 7 can receive sizable contributions from the new 
charged scalars 77^^, |[^IT7ll4^1 - l5^ . The partial decay widths for these processes can be expressed 

in terms of the following ratios: 


(Tii (pp —>■/i —)• 77 ) (Tii(pp—)-/i) BRn((i—)- 77 ) 

o'sm(pp crsM(pp h ) BRsm(/i 77) ’ 


and similarly for Rz'y- The current signal strength in the channel pp —)• /i —)• 77 is 1.17±0.27 at ATLAS 
ll54l and 1 . 141 q 23 CMS ll55l . Both ATLAS and CMS are consistent with SM at Icr level, but still have 








a relatively large uncertainty. Detailed analysis shows that the signal rate is totally determined by the 
sign of A4 and A4 + A 5 / 2 . is enhanced (suppressed) for negative (positive) A4 with a small variation 
caused by A 5 Il50ll53l as A 5 is tightly constrained by the electroweak precision data. Ref. |I52]| shows that for 
a given enhancement in the decay rate, an upper bound can be set on the type II seesaw scale; for instance, 
for an enhancement of 10%, the upper limit on the seesaw scale is about 450 GeV, which is completely 
within the reach of the 14 TeV LHC. With this conclusion, dedicated searches at the LHC for the decays 
and —)• H^W^* with Mu±± < 450 GeV would be strongly motivated. 

The current collider limits on the type II seesaw are mainly drawn from the doubly-charged scalar 
search through its like-sign dilepton signature, —)■ At LHC 7 TeV, L = 4.7 fb“^ and assuming 

BR(H=^=^ ^ ~ 100%, ATLAS EH excluded Mh±± below 409, 398, 375 GeV at 95% C.L. in the 

e^e^, channel respectively. With the same assumption, CMS |[57l set a lower bound ranging 

from 204 GeV to 459 GeV. Recently, ATLAS has updated its limits to LHC 8 TeV, L = 20.3 fb~^, and 
pushed the most stringent lower limit up to 550 GeV 15^ . We should emphasize that these bounds do not 
apply to the nondegenerate case especially in the negative scenario, where two other decay modes of 
the like-sign diboson decay —)• and the cascade decay —)■ can dominate 

over the like-sign dilepton decays in a large portion of parameter space. To our knowledge, there are no 
direct experimental limits for these two channels so far. For the like-sign diboson decay. Refs. 11591 - 1611 
concluded that the lower bound on Mfj±± can be derived from the like-sign dilepton ATLAS limit, which 
is 60 GeV for LHC 7 TeV, L = 4.7 fb“^, and extends to 84 GeV for LHC 8 TeV, L = 20.3 fb“^. We 
will show later that the constraints from the cascade decay may be comparable with the one obtained from 
the like-sign diboson signature. As was also shown in Fig. 5 of EH, the mass limit on depends 
significantly on the branching ratios assumed for the like-sign dilepton decays, and degrades quickly with 
a decreasing branching ratio. For instance, the benchmark point chosen for our simulation in the positive 
scenario Mjj±± = 400 GeV is located in the safe region, as long as BR(i7=*==*= ^ < 0.5. From 

Table shown in Sec. IV A we see that a potential danger may occur in the channel 




e e 


which is almost at the edge of the exclusion limit. In the negative scenario, we can choose a much lighter 


as our benchmark point, see Eq. (301, because its cascade decays can dominate overwhelmingly over 
the dilepton modes in a large portion of parameter space; for instance, it can be lighter than 150 GeV if 
—)■ < 0.02. In this case, the two triplet-dominating neutral scalars H^, can be nearly 

degenerate with the SM-like Higgs h, and could result in interesting phenomena. To put it in short, taking 
into account the rich decay modes in the nondegenerate case, our benchmark points are compatible with the 
most recent ATLAS constraints on type II seesaw. 


The previous direct searches at LEP also put constraints on new scalars. Eor the neutral ones, EEP set a 
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lower bound m^o > 80 — 90 GeV for general models with an extended scalar sector ll62ll . For the charged 
scalars, searches were performed for pair production; for instance, in minimal supersymmetric standard 
model (MSSM) or in the more general type-II two-Higgs-doublet model (2HDM), the combined LEP data 
yielded mu± > 80 GeV or 72.5 GeV 1631. Although these lower bounds on masses are generally model 
dependent, they are respected in our later numerical analysis. 


III. DECAY PROPERTIES OF THE SCALARS WITH NONDEGENERATE MASSES 

The decay properties of the new scalars have been studied in the previous literature ifTdl |33l |34l HH ; 
see, in particular. Refs ll^ |34l for the relevant formulae of the decay widths. Since these properties will 
be important in devising signal channels at colliders, we summarize them in this section and present a few 
figures that would help the reader understand our later analysis more readily. We are aware that some similar 
figures have been ploffed in fhe above liferafure albeif for differenf paramefers fhan ours. 

As we menfioned earlier, fhe cascade decay —)■ or —)■ is possible in 

fhe nondegenerafe case, which can change significanfly fhe decay paffems even for a small mass spliffing 
AM. These decays are classified info fhree cafegories: lepfonic, gauge boson, and cascade decays, and fheir 
branching rafios are respectively mosf sensifive fo fhe Yukawa coupling Y, friplef vev V/\, and mass spliffing 
AM. From fhe decay phase diagram presenfed in Ref. Il64l . one knows fhaf lepfonic decays dominafe for 
bofh AM and va being small, gauge boson decays dominafe in fhe region of small AM and large va, and 
cascade decays are dominanf for a large AM and a moderate va- In parficular, for va around 10“^ GeV, 
cascade decays dominate in a large porfion of fhe AM paramefer space even if AM is as low as 2 GeV. 


A. Positive Scenario 

In fhe positive scenario, Mjj±± < Mfj± < Mho ^^o, fhe doubly charged scalar is fhe lighfesf and 
fhus shares decay modes similar fo fhe degenerafe case, i.e., fhe like-sign dilepfon decay if if 

{i = e, /i, r) and diboson decay —)• Their decay amplifudes are proporfional fo fhe Yukawa 

coupling Y and fhe friplef vev va, respecfively. Figure shows fheir branching rafios as a funcfion of va 
af Mff++ = 400 GeV on fhe leff panel, and as a funcfion of Mfj++ af va = 10““^ GeV on fhe righf. 
One can see fhaf fhe fwo decays are comparable around va ~ 10“^ GeV and M^++ = 400 GeV. For 
a given Mh++, fhe dilepfon decay dominates af va < 10“®GeV, while fhe diboson decay lakes over af 
va > 10“'^ GeV; conversely, given va, the dilepfon (diboson) decay dominates for a lighter (heavier) 
A plof similar fo fhe leff panel of Fig. [^can be found in Refs. ll^[3^ af a lower mass. 
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VA(GeV) 



MH++(GeV) 


FIG. 2. Decay branching ratios of as a function of va (left panel) and M^++ (right) in the positive scenario. 

The heaviest scalars ijj^ = H^, follow the cascade decay chain —)■ H^W^* with 

The branching ratios are plotted in Fig. |^as a function of va at AM = 10, 30 GeV and 
Mjj++ = 400 GeV. We notice some interesting features. First, for a given AM, the cascade decays 
of H'^, and share a similar region and shape in the plots. Second, for AM = 10 GeV, the 
cascade decays dominate in the region 10“^ GeV < fA < 10“^ GeV. The leptonic decays dominate for 
VA < 10“® GeV, which is about three orders of magnitude lower than in the degenerate case. On the other 
hand, gauge boson decays dominate when va > 10“^ GeV, which is three orders of magnitude higher than 
in the degenerate case. For AM = 30 GeV, the domination region of cascade decays is further extended to 
10“® GeV < VA < 10“^ GeV, making the other two modes less important in a larger parameter region. 

B. Negative Scenario 

We choose the following benchmark masses to illustrate our results in the negative scenario: 

AM = 10 GeV : M^++ = 149 GeV, Mh+ = 140 GeV, Mho^a^ = 130 GeV; 

AAF = 30 GeV : Mj^++ = 185 GeV, = 160 GeV, ~ 130 GeV, (16) 

where Mfj++ is worked out from M^o and AM by Eq. Q. They will also be used in our later collider 
simulations in SecIVl 

The decay branching ratios of H^~^, H^, A^ are shown in Fig. as a function of va- The 

same plots were presented previously in Refs, ll^ for either smaller or larger masses. The decay 
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FIG. 3. Decay branching ratios of 7F+ (upper panel), (middle), and (lower) as a function of va in the positive 
scenario, with Mh++ = 400 GeV and AM = 10 GeV (left panel) and AM = 30 GeV (right). 



















12 


properties of have changed significantly from the degenerate case. The cascade decay 

is dominant in a large region of parameters. For instance, —)• dominates in the 

range GeV < ua < 1 GeV (10"^ GeV < r;A < 1 GeV) for AM = 10 GeV (AM = 30 GeV), 
while the like-sign dilepton and diboson decays are heavily suppressed compared to the degenerate case (see 
Fig.[^. Actually, the like-sign dilepton decays can be safely ignored when va > 10 ® GeV (10 ^ GeV) 
for AM = 10 GeV (30 GeV), thus in the majority of parameter space that we are interested in, the current 
LHC bound on can be easily avoided. The dominant regions for —)■ , A^W~^* are similar 

to those of 

For the neutral scalars, the relevant decay modes are —)• bb, r+r", cc, W^W~, ZZ, hh and 

> bb, T+r”, cc, Zh. (The branching ratios for —)■ hh and A^ —)• Zh are too small to show 

properly in Fig.|^) The invisible decay H^/A^ —)■ uv dominates in the small va region, while the hadronic 
decay /AP —)> bb dominates for a larger va- When va > lO”"^ GeV, — )• can reach 


about 10%, which is a useful channel for probing neutral Higgs scalars as we will discuss in Sec. V C 


IV. LHC SIGNATURES IN THE POSITIVE SCENARIO 

The main goal of this section and the next is to investigate the LHC signatures of type-II seesaw with a 
nondegenerate spectrum. Since there are many possible final states resulting from cascade decays of scalars 
and the decay patterns are also different for the positive and negative scenarios, we treat them separately. 
Before going to the details, we briefly summarize our simulation procedures. We implement the model in the 
Mathematica package FeynRules ||65l, whose output UFO ll^ model file is taken by MadGraphS ETl 
to generate parton level events for relevant physical processes. Those events then pass through Pythia6 
ll^ to include the initial- and final-state radiation, fragmentation, and hadronization. Delphes3 |[69l is 
then used for detector simulation and MadAnalysisS f/Oll for analysis. In our simulation, we choose 
to work with the CTEQ6L1 parton distribution function (PDF) iTll . The tree-level total production cross 
sections of scalars at 14 TeV LHC (LHC 14) have been plotted previously as a function of scalar masses for 

channels in Ref. lfT4l and for the H^A^, channels in 

Refs. l|33l[33|. The interested reader should refer to those papers for details. 

Some of the signal channels to be studied here were considered in previous papers ESi 13X113311 . but most 
of those analyses were based on theoretical estimation or parton level simulation and thus cannot be directly 
compared with realistic experimental data. In contrast, we simulate the signal channels for both positive 

' In this work, we consider only the tree-level contributions. The QCD correction to H~~ pair production was computed 
in m, with a A-factor of about 1.25, while the contribution from real photon annihilation tends to increase the production by 
10%, resulting in an overall if-factor of 1.35 Il4l . The H~ associated production in principle gives a similar A-factor 
~ 1.25. 
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FIG. 4. Decay branching ratios of (upper panel), (middle), (lower right) and (lower left) as a 
function of va in the negative scenario, with M^o = 130 GeV and AM = 10 GeV (left panel) and AM = 30 GeV 
(right). being the lightest, their branching ratios are independent of AM. 
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and negative scenarios at the detector level and design specific cut criteria for each channel, and this will 
be more powerful in signal prediction. In the following subsections we present our analysis for each signal 
channel in the positive scenario, and devote the next section for the negative scenario. To facilitate lepton 
and charge identification, the leptons here refer only to electrons and muons. 

A. Signals for doubly-charged scalars 

In the positive scenario, the doubly-charged scalars are the lightest and decay directly into SM 
particles. Therefore, their signatures are essentially analogous to the degenerate case, the most promising 
signal being still the four-lepton channel: 

pp H+^H— £+£+ -h rr. (17) 

The branching ratios of the dilepton decays depend on the lepton flavor and neutrino mass hierarchy. They 
can be easily worked out ifldll . and the numbers are shown in TableWe observe the following relations, 

e^e^) > > BR{H^^ t^t^) for IH, 

BR{H^^ ^ T^T^) >BR{H^^ ^ p^p^):$>BR{H^^ ^ e^e^) for NH, (18) 

which arise as a consequence of the mass hierarchy and mixing pattern ifTSlfT^IT^ . The main irreducible 
background for this signal channel is 

zz ^ i+r +£+r, (19) 

and reducible backgrounds are mainly from tt, Zbb, Ztt. We simulate the channel at 14 TeV LHC with an 
integrated luminosity of 300 fb~^ (LHC14@300) and assume va = 10“^ GeV, AM = 30 GeV. As we 
mentioned in Sec.|T^ the latest ATLAS bound on Mu±± can be relaxed to ~ 400 GeV for a nondegenerate 
spectrum, we thus choose the four values of M^±± = 400, 500, 600, 700 GeV for both NH and IH. To 

perform more realistic simulation, we use the same selection criteria as in the ATLAS paper Il56ll . 



ee 

efi 

er 

pp 

pr 

rr 

IH 

48.7 

0.498 

0.602 

14.6 

24.5 

11.1 

NH 

0.793 

4.23 

0.177 

27.6 

30.7 

36.5 


TABLE 1. Dileptonic branching ratios (in percentage) of to different flavors. 

The distributions in the like-sign dilepton invariant mass Mi±£± for the signal and background are plot¬ 
ted in Fig. 1^ for both NH (left panel) and IH (right) cases. We see that the distributions of signal and 
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FIG. 5. Reconstruction of via dilepton invariant mass M(±i± at LHC14@300 and for Mfj±± = 

400, 500, 600, 700 GeV. 
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FIG. 6. Same as Fig.|^but upon imposing the cut ( |20| i. 

background are well separated: compared with the signal, the background is mainly located in the small 
Mi±i± region. Moreover, the IH case has larger signal events, because have larger branching ratios 
to the electrons and muons (see Eq. ( [TS] )) which are counted as leptons in simulation. To further purify the 
signal, we impose the cut 


Mi±i± > 300 GeV. (20) 

In Fig. we show the distributions of Mi±£± in the range 300 — 800 GeV after imposing the cut. As 
an illustration of details. Table [T^ shows the evolution of survival numbers of events, statistical signifi¬ 
cance S/\/S + B, and signal to background ratio S/B upon imposing the cuts step by step for Mj£±± = 
400 GeV. Here the first three cuts are just the ATLAS ones ||56ll . while the final one is Eq. ( [20] ). One can 
see fhaf fhe backgrounds now become negligible. As a consequence, af LHC14@300, one has fhe pofenfial 
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to probe the doubly-charged scalars with a mass up to ~ 600 GeV for NH and ~ 700 GeV for IH. 


cuts 

signal 2 f± 2 f=F 

ZZ 

S/^/S + B 

S/B 


IH 

NH 


IH 

NH 

IH 

NH 

no cuts 

468 

136 

29276 

2.71 

0.798 

0.016 

0.00467 

> 20GeV,p|, > 25GeV, Ai?« > 0.4 

446 

131 

15921 

3.49 

1.03 

0.0280 

0.00820 

Mi±^± > 15GeV 

446 

131 

15919 

3.49 

1.03 

0.0280 

0.00820 

70 < M^±^± < llOGeV (Zveto) 

446 

131 

10200 

4.33 

1.28 

0.0438 

0.0128 

Mt,±(± > 300GeV 

194 

56.8 

3.19 

13.8 

7.33 

60.8 

17.8 


TABLE 11. Survival numbers of four-lepton signal from with Mh±± = 400 GeV and its main background 

ZZ, statistical significance S/\/S + B, and signal to background ratio S/B at LHC14@300. 


B. Signals for singly-charged scalars 

To illustrate our analysis for the production and detection of singly-charged scalars, we choose M}j±± = 
400 GeV and AM = 30 GeV so that = 458 GeV. In the positive scenario, cascade decays 

into via the radiation of an off-shell . Here we consider the five-lepton signal coming from 

associated production 0 

pp -h -h t^T-, (21) 

where the off-shell W^* decays leptonically. The main SM background is 

ZZW^ tt- + .(+£- + (22) 

The simulation is performed at LHC14@300. We apply the same cuts as for the four-lepton signal. The 
numbers of events before and after cut selections for both signal and backgrounds are listed in Table H 
Since BR(Z —£~^£~) is only about 6%, the background ZZW^ is already much smaller than the signal at 
the pre-selection level. After imposing the cuts, we have 27.2 signal events for NH case, with a statistical 
significance of 4.73. So for the NH case, a 430 GeV singly-charged scalar can be actually discovered 
if we consider additional contributions from H^A^, and production. The IH case is even 

better. According to our simulation, 104 signal events pass the selection cuts, with a statistical significance 
of 9.92. Thus a 5cr significance discovery requires L = 76 fb~^ at LHC14. 


^ The five-lepton signal can also originate from H^H^, and associated production 1281 . which can potentially 

enhance the signal strength. We do not consider these additional channels for simplicity. 
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Pre-selection 

Post-selection 

S/^/S + B 

S/B 

NH 

40.5 

27.2 

4.73 

4.63 

IH 

154 

104 

9.92 

17.7 

ZZW^ 

9.00 

5.87 

— 

— 


TABLE III. Numbers of events and statistical significance S/\/S + B, S/B before and after imposing cuts for 
Mh±± = 400 GeV, Mh± = 430 GeV at LHC14@300. 


The final states originating from decays include missing particles, and therefore we cannot fully 
reconstruct its invariant mass. Since the leptons coming directly from decays are generically much 
more energetic than those from off-shell W decays, we can employ this feature to design a kinematical vari¬ 
able. Consider the cluster formed by a like-sign dilepton and the least energetic lepton of opposite 

charge ) and define the cluster transverse mass, 

^ +*)' - +*)'. ( 23 ) 

where are respectively the transverse momentum and invariant mass of the 

cluster. [Similar notations will be used below.] The peak structure of this variable then indicates the mass 
of as can be clearly seen in its distribution in Fig.j^for both IH and NH cases. 



FIG. 7. Reconstruction of via for Mh± = 430 GeV at LHC14@300 
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C. Signals for neutral scalars 

Now we consider the heaviest particles in the positive scenario, the neutral scalars , AP. The previ¬ 
ous literature has discussed possible like-sign four-lepton signals from associated production involving the 
neutral scalars and H^A^) which cascade decay into like-sign doubly-charged scalar pair 

plus far off-shell W*, with being detected in their dilepton channel 12814301 . However, to 
gain sizable signal events, the mass splitting has to be limited in AM ~ 1 — 4 GeV OOl . Such a splitting 
may be too small to be directly detectable at a hadron collider, and this signal channel may only be consid¬ 
ered as indirect evidence for neutral scalars. In this work, we consider a more promising channel, i.e., the 
LNV signal: 


pp (24) 

where -0° = A^, and denotes leptons from off-shell W decays, and the two jets in the final state 

are too soft to be detected. 

We emphasize that the above signal channel has distinct features compared with previous studies. For 
instance. Refs, ll^l^ focused on another LNV channel: 

pp ^ (25) 


They also noted that due to significant interference between Hq and for ~ M^o, the usual zero 
width approximation does not apply. The two channels can be labeled by the final leptons from the 
decay plus off-shell VF*s, with the cross sections being 


a{2i^2t^ + 2 W*^W*^) = a{pp , 

a{Ai^ + 3W*^) = a{pp , (26) 


where BRcas. = BR(M=*= —)■ H^^W*^)BR{H^^ —)> ifif)- One sees that the only difference between 
the two channels is the cross section a{pp —?■ H^W^*) versus a{pp ——)• 

Following the generalized narrow width approximation introduced in Refs, f/dll . we have calcu¬ 
lated the above cross sections including the interference effect, with the result: 


a{2l^2t^ + 2W*^W*^) 
a{U^ + 3VF*^) 


a{pp —^ 
a{pp —)• 


BR^o “h BR^o “h 2 
BR^o BRyi^o — 2 


BR^o+BR^q V 
BR^o BRy^o J 


BR^o +BR^o Y 
BRji^o -f BRy^o J 


BR 


2 

cas.! 


BR 


2 

cas. 1 


(27) 
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where BR^o/^o = —)• For rigorously degenerate masses and widths, M^o = M^o 

and BRjyo = BR^o, Eq. ( [T/] ) reduces to 

a{2l^2t^ + 2W*^W*^) = 2a{pp ■ 2BR^o/^o • BR^^^., 

a{Al^ + 3PF*^) = 0. (28) 


Therefore, upon taking into account the interference effect between the intermediate Hq and states. 


our signal in Eq. (24i is enhanced by approximately a factor of two while the signal in Eq. (25 I tends to 
disappear. We have confirmed this effect by Madgraph simulation, and the vanishing result for the process 
( [25] ) is also consistent with the statements in Refs. li30ll3Tl . 

Our signal violates lepton number by two units and thus has no irreducible SM background, but its cross 
section is also relatively small. Eor EHC14@300, there are only 11.4 signal events for NH and 38.0 for 
IH. We therefore do not apply any cuts to it. We combine the like-sign dilepton and the two least energetic 
leptons of opposite charge into a cluster and define fhe clusfer fransverse mass: 








+Xt^ 








(29) 


Ifs disfribufion ploffed in Eig.j^can be employed fo parfially reconsfrucf fhe neufral scalars af M^o = 
458 GeV. 



FIG. 8. Reconstruction of 72°, A° via for Mho/ao = 458 GeV at LHC14@300. 


V. LHC SIGNATURES IN THE NEGATIVE SCENARIO 


In the negative scenario, Mfj±± > Mu± > M^o ^o, the more charged scalars can dominantly cascade 
decay into less charged ones with the radiation of an off-shell boson in a significant portion of parameter 
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space. This makes collider search strategies of the scalars very different from the degenerate case, and a 
systematic study is still lacking in the literature. In this section we shall investigate the issue and when 
possible compare our results with those in the literature. For illustration, we assume, unless otherwise 
stated, Mjjo = 130 GeV, AM = 30 GeV, so that our benchmark point for the scalar masses is 

Mjj±± = 185 GeV, Mfj± = 160 GeV, Mjjo y^o = 130 GeV. (30) 

A. Signals for doubly-charged scalars 

The doubly-charged scalars are the heaviest in the negative scenario, and some signal channels for 
their cascade decays have been briefly discussed in Refs. ll^l64l . Considering that the cross section for 
the pair production {(Tfj++H — ~ 98.5 fb at LHC14) is smaller than that for the associated production 
{aH±±HT ~ 223 fb), we concentrate on the latter with bb^ 

pp i^i^^Tbb + bbjj. (31) 

In this signal, the like-sign W^* pair decays leptonically while the other oppositely charged W^* decays 
hadronically. Since the jets from the latter are very soft and difficult to detect at LHC, the visible final sfafes 
appear as i^i^^Tbb -I- bb. The like-sign dilepfon will make fhe channel suffer from less backgrounds. 

We sfarf wifh some basic cufs: 

p^ > 10 GeV, \rji\ < 2.5, 

pip > 20 GeV, \r]j\ < 2.5, > 30GeV, 

ARu > 0.4, ARjj > 0.4, ARje > 0.4. (32) 

The leading irreducible background is 

ttW^ W+bW-bW^ i^i^^rbbjj, (33) 

while ofher backgrounds like W^W^jjjj, tb are much smaller |[T4l . To keep fhe signal rale as much as 

possible, we do nol apply 6-fagging; insfead, we impose fhe evenl selecfion = 2 lo pick up events with 

exactly a like-sign lepton pair. Furthermore, since both leptons are from off-shell W s, we apply a veto cut 
on the largest transverse momentum of the lepton pair, pp < 30 GeV. 

Figure displays the distributions of particle separations AR^^ ji after imposing above cuts on both 
signal and tiW^ background. It is evident that the leptons and jets from the ttW^ background are more 


^ This signal can also be produced in the H pair production. A rough estimation based on the ratio of cross sections 
shows that it enhances the signal events by a factor 1.5. 
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isolated than those from the signal. We can thus further purify the signal by demanding 

iS.Ri±i± < 2 . 0 , < 2 . 0 . 


( 34 ) 


We list in Table IV the survival numbers of events, statistical significance, and signal to background ratio 
upon imposing the cuts. We see that all the cuts chosen here, especially those on and Ai?, are efficient 
enough. The statistical significance can reach 4, and we have about 16.5 signal events at LHC14@300. 




FIG. 9. Distributions of particle separations (left panel) and (right) after imposing the cuts Ni± = 2, 

< 30 GeV, and Eq. at LHC14. 


Cuts 


ttw^ 

S/^/S + B 

S/B 

Basic Cuts 

549 

3682 

8.45 

0.1463 

N,±=2 

70.4 

1588 

1.73 

0.0444 

< 30 GeV 

50.0 

63.8 

4.69 

0.785 

/S.Ri±i± < 2.0 

40.2 

23.6 

5.04 

1.71 

ARji < 2.0 

16.5 

1.35 

3.91 

12.2 


TABLE IV. Survival numbers of signal and background, statistical significance S/^/S + B, and signal 

to background ratio S/B upon imposing each cut sequentially at LHC14@300. 


All new scalars appear in this signal channel. The neutral ones , ^4*^ decay with no missing particles 
and can be fully reconstructed by using the two 6-jet invariant mass Mjj. The charged scalars 
decay with missing energy and can be partially reconstructed with the help of the cluster transverse mass: 


= 


PtJJU + ^jjU - {pTJjU +Xt 


(35) 


The distributions of these two variables are plotted in Fig. 10 One can see the peaks at Mjjo = 130 GeV 


and Mjj±± = 185 GeV respectively, although the numbers of events are limited. 
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FIG. 10. Distributions in Mjj and in signal at LHC 14 @300. 


Here we would like to comment on some similarities and differences between cascade decay signal from 
in our study and like-sign diboson signal from in previous studies |[T3llT4ll . First, if 

we do not use 6-tagging, the signals are the same for both cases. But in practice, the jets are 

quite different. In the cascade decay, all four jets are 6-jets, which come from decays of the neutral scalars 
A^, while in the diboson case, all four jets are light ones from decays of like-sign iy=*=*s. The invariant 
mass Mjj can be used to distinguish between these two signals, since M^ois usually heavier than W^. 
Also considered in Ref. |[T4l are the contributions from gauge decays of which could produce 

signal So according to the number of 6-jets, the original sources of the jets system can be 

clearly distinguished from each other once 6-tagging is applied. 

Next, we look more closely into the like-sign dilepton signal. We note that the like-sign dilepton from 
cascade decays of in our study is similar to that from leptonic decays of like-sign di-VF in the pair 
production of but with a lower mass llT3l : 

cascade decay: pp —> > H^W^* -t- —)■ + jj{jj), (36) 

diboson decay: pp W^W^* -f (37) 


Again, {jj) in Eq.(361 denotes undetectable soft jets. For a comparative study, we use the masses as shown 
in Eq. for our cascade decay channel, while in the diboson channel a lower mass Mjj±± = 90 GeV 
(case A) or Mh±± = 150 GeV (case B) is assumed as in Ref. llT3ll . The distributions of some kinematical 


variables with no cuts are depicted in Eig. 11 


We note some features in the distributions. Eirst, the leading lepton transverse momentum p^ is quite 
different in the three cases. In the diboson decay, a heavier Mfj±± usually implies a larger pf^, while p^ 
from cascade decay is generally much softer. Eor the next-to-leading lepton transverse momentum p^, the 
cascade decay and diboson decay in case A are similar. Second, the distributions of are similar in 
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FIG. 11. Distributions of ARuj^, and Mu in cascade decay H^W^* and like-sign diboson 

decay —>■ W^W^* in the signal event at LHC14. 


all three cases, but the distributions of can be very different. This difference arises mainly because 
in cascade decay the dilepton and two of four jets are from one single and the other two jets from 
t/^o> resulting in a ARji with two distinct peaks, while the dilepton and jets in diboson decay come from 
different and are thus more isolated with a single peak in ARj£. Third, the distributions in missing 

transverse momentum are almost the same in cascade decay and diboson decay (case A). Finally, the 
distributions of dilepton invariant mass Mu are also distinguishable, with a peak around 20, 30, 70 GeV 
in cascade decay, case A and B, respectively. To summarize in a word, the distributions of p^, ARji, and 
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Mu can be used to distinguish among these three cases of signals. 

B. Signals for singly-charged scalars 

The LHC physics of the singly-charged scalars in type-II seesaw has been studied in many previous 
papers but only for the degenerate case, see, e.g.. Refs. llT4l ITtI [TSll . In this subsection, we carry out a 
detailed study for the nondegenerate case in the negative scenario. As one can see in Fig. the cascade 
decay —)• = H^, A^) dominates in a wide interval of the triplet vev, 10“® GeV < < 

1 GeV. For va. ^ 10“^ GeV, decays to vv, more interesting is the case va ^ 10“^ GeV, when 
decays mainly to bb, r+r”, and to bb, T+r“, W~^W~. We therefore consider first the following 
channel 0 

pp ^ + bbbb. (38) 

Since the charged lepton and neutrino come from an off-shell W, their energies are relatively smaller than 
those from an on-shell W in the SM background. This feature can be utilized to improve the signal to 
background ratio efficiently. 

The main SM backgrounds we considered are: 

ti ^W^bW-b^i^^T + bbjj, (39) 

tibb ^W+bW-bbb + bbbb + j j /l^ , (40) 

W^bbbb^ + bbbb, (41) 

where the first two are reducible and the last one is irreducible background. For the tt background, we only 
have to consider the semi-leptonic decay channel with the two light jets mis-tagged as 6-jets. For ttbb, both 
semi-leptonic and di-leptonic decays of tt are involved. In the di-leptonic decay channel, we require that 
one of the charged leptons escape detection. The identification efficiencies for e and p are taken to be 0.9 
and 0.95 respectively, and L = 30 fb~^ is assumed at 14 TeV LHC (LHC14@30). 

To isolate the signal, we employ the following basic cuts: 

p^ > 10 GeV, 1%! < 2.5, 
p^P > 20 GeV, \rij\ < 2.5, 

ARei > 0.4, ARji > 0.4, > 0.4. (42) 

We notice that the same l^Mtrbhbb signal could also appear in MSSM via the /AP associated production 1751 : 

pp tb/ib + bb^ W^bbbb l^^rblbb, 

where in this case Mu± is heavier than the top quark. Besides different thresholds of H^, the decay modes of can also 
he used to distinguish between the two scenarios. One distinct feature is that the W bosons in our cascade decays are always 
off-shell while those in MSSM from top decays are on-shell. Another one concerns the separation of the two 6-jets in 
decays. In cascade decays, the 6-jets originate from /A^ decays, while in MSSM one comes directly from and the other 
from tops, resulting in generally more separated two 6-jets than in cascade decays. 
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Cuts 


tt 

ttbb 

W^bbbb 


S/B 

S/VS + B 

Basic Cuts 

567 

2039202 

162410 

1295 

0.373 

0.000256 

II 

o' 

II 

II 

22.2 

317 

218 

3.45 

0.940 

0.0412 

p^j, < 30 GeV 

18.4 

69.6 

53.7 

0.488 

1.54 

0.148 

Xt < 30 GeV 

10.6 

5.53 

3.80 

0.0916 

2.37 

1.13 

M|. < 30 GeV 

9.00 

1.63 

1.98 

0.0305 

2.53 

2.47 

< 2.0 

8.94 

1.30 

1.32 

0.0305 

2.63 

3.37 


TABLE V. Survival numbers of events, statistical significance S/^/S + B and signal to background ratio S/B after 
imposing each cut sequentially at LHC14@30. 

We found that for Mu± — M^o = 30 GeV, the efficiency of basic cuts for a signal lepton is about 0.77. To 
further purify the signal, we first apply the event selection: 

Nb = 4, Nj = 0, Ne = 1, (43) 

which means that we only choose the events which exactly contain four 6-jets and one charged lepton. We 
found that better 5-tagging efficiency could potentially suppress the leading background coming from the tt 
channel. The next to leading background is ttbb channel, while the events from W^bbbb are already smaller 
than the signal. 

To suppress further the background, we employ several kinematical cuts that are largely designed on 
the off-shell nature of the W in the signal versus the on-shell nature in the background. The various dis¬ 
tributions are displayed in Figj^ after imposing the basic cuts. The transverse momentum of leptons 
in all backgrounds has a peak at about 30 ~ 40 GeV, while in the signal is much softer. The distri¬ 
butions of the missing transverse energy share a similar feature. The off-shell or on-shell origin of 
leptons and neutrinos can be directly probed in the lepton transverse mass distribution. Since the mass 
splitting between and H^/A^ is 30 GeV, most signal events have < 30 GeV. In contrast, the 
background events have a clear edge around ~ 80 GeV and only a small number of them falls in the 
region < 30 GeV. Finally, the distribution in ARb£ has two peaks for the signal, one around 0.9 and 
the other around 3. The smaller one most likely corresponds to the case when both particles are from the 
same cascade decay —)• —)■ i^ubb, while the larger one has the 5-jet more likely from the decay 

of the directly produced The separation in the background is usually large and has a peak at 3. These 
differences can be used to distinguish between the signal and backgrounds, through the following cuts: 


p^T < 30 GeV, Xt < 30 GeV, < 30 GeV, ARbi < 2.0. 


(44) 
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PT(GeV) (GeV) 




FIG. 12. Distribution of pfp,^T, Mip, and ARte after imposing basic cuts for signal i^^Tbbbb and its backgrounds. 


Table [v| shows the survival numbers of events for the signal i^^Tbbbb and backgrounds after the basic 
cuts and further sequential cuts, together with S/^/S + B and S/B. The cuts are efficient enough in 
preserving the signal while suppressing the background. At LHC14@30, about 9 signal events survive 
with a statistical significance of 2.63 and a signal to background ratio of 3.37. Then for a 5a reach, the 


required luminosity is 109 fb“^. As shown in Fig. 13 the mass of H^/A^ can be fully reconstructed using 
the invariant mass of the two f)-jets in the production. Since and A^ are degenerate to good 

precision, we see effectively one peak at 130 GeV in the distribution of Mi,b- To reconstruct the mass of 
we still use the cluster transverse mass: 


= 


PT,bbi + 


Mlt+Jiry 


PT,bbe +^t) , 


(45) 


whose distribution is shown on the right panel of Fig. 13 


For the associated production, we consider the second channel: 

pp ^ + bbr+T-. 


(46) 


This process has a much smaller rate than the previous one in Eq. ( [38] ), but its background is also tiny. We 
thus use LHC14@300 as an illustration. The ft/r-tagging efficiency and mis-tagging rate are assumed to be 
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FIG. 13. Distributions of and for i^^xbbbb signal channel at LHC14@30. 
the same as in our previous analysis. The main backgrounds are: 


ttT~^T 

W~^bW bT~^T + bbT~^T +jj/t^, 

(47) 

tiw^ 

W+bW-bW^ + bbr+T-, 

(48) 

W^bbT^T~ 

"—>■ + bbT'^T~ . 

(49) 


For the tiT~^T~ background, we consider both semi-leptonic and di-leptonic decays of tt. One of the 
charged leptons in the di-leptonic channel is required to escape detection, while the additional two light jets 
in the semi-leptonic channel receive no further constraints. The two irreducible backgrounds ttW^ and 
W^bbT+T- are also taken into account in our study. As in our previous analysis, we start from basic cuts. 


We apply the same basic cuts as in Eqn. (42 1 for the signal channel i^^xbbbb. The evolution of the event 


numbers and S/\^S + B and S/B upon each cut is summarized in Table VI At the level of the basic cuts, 
the main background comes from tiT~^T~, while the two irreducible backgrounds are already smaller than 
the signal. 


Cuts 


Ut^t 

-H 

■40 

W^bbr'^T 


S/B 

S/^S + B 

Basic Cuts 

1905 

9935 

356 

305 

16.28 

0.171 

Vf, = 2,A, =2,A^ = 1 

62.8 

76.1 

8.22 

3.36 

5.12 

0.717 

110GeV< Mfeb <150 GeV 

39.7 

14.3 

1.53 

0.300 

5.31 

2.46 

A.Ri)h < 2.0 

27.7 

5.85 

0.450 

0.0900 

4.74 

4.30 


TABLE VI. Survival numbers of signal i^^Tbbr'^T and its backgrounds, statistical significance S/^/S + B, and 
signal to backgrounds ratio S/B after imposing each cut sequentially at LHC14@300. 
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We then apply the event cuts: 


Nk = 2,Nr = 2,N^ = l, (50) 

to select the signal . The numbers of signal and background events are now of the same order. 

In contrast to our previous signal channel l^^xbbbb, we see that requiring two r-jets from four jets final 
states are very useful to remove the backgrounds. When comparing to the signal channel in the 

production, the additional charged lepton is also very helpful to suppress the backgrounds. Except 
for a smaller production rate, the signal £^^TbbT~^T~ is the cleanest among what we have investigated. 




FIG. 14. Distributions of Mi,b and Mrr for i^^Tbbr'^T 


In Fig. 14 we show the distributions of invariant masses M^b and M^t after implementing the cuts 


(50l. There is a clear peak at 130 GeV for the signal in the Mbb distribution, which can be employed to 
reconstruct the neutral scalars But the M^t distribution is not distinctive enough. There are two 

reasons for this difference. First, the r pair in the background mainly comes from the on-shell Z decay, 
whose mass is close to that of the neutral scalars; second, the peak of M^t is shifted to the low mass region 
compared with the original value since we do not renormalize the momentum of the r jets. To reconstruct 
the charged scalar, we need implement more cuts. We first narrow down to a window of Mbb- 


no GeV < Mbb < 150 GeV. 


(51) 


Next we concentrate on the off-shell W boson again. In Fig. 15 the distributions in p^, and M^ are 
plotted. They are changed significantly from those of the previous signal i^^Tbbbb. The leptonic decays 
of the T lepton contribute to the final sfafe lepfons 1. These fake lepfons have a much harder fransverse 
momenfum and behave jusf like fhe backgrounds when > 20 GeV. The missing fransverse momen- 
fum is also much larger, since additional (anfi-)neufrinos carry away parf of fhe energy in fhe hadronic 
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decays of the r leptons. This also alters the distribution of the lepton’s transverse mass . Even worse is 
that the distribution of happens to be similar to that for the largest ttT~^T~ background. Therefore, the 
cuts applied on the off-shell W^* would not be very effective. 






FIG. 15. Distributions of and Ai?;,;, in signal at LHC14. 

Also shown in Fig. [T^is the distribution in separation AR^b- The two 6-jets in the signal come from the 
decay of one heavy neutral scalar while those in the background originate from the decays of tt respectively. 
This results in a smaller ARbb for the signal than for the background. So we apply one more cut: 

ARbb < 2.0. (52) 

After all the cuts, we still have about 28 signal events and 6 background events, with a statistical significance 
around 4.7 and a signal to background ratio 4.3. This signature is therefore also promising at LHC14@300. 

Theoretically, the mass of could be reconstructed using some transverse variables of the bM or rri 
system. But since has additional contributions from hadronic r decays, the bb£ system is not useful 
here. We therefore employ the cluster transverse mass for the tt£ system: 

, ( 53 ) 

whose distribution shown in Fig. [^displays a clear peak albeit with a small number of events. 
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FIG. 16. Disti-ibutions of for i^^Tbbr+T- signal channel at LHC14@300. 


Finally, we would like to discuss briefly some other interesting signals in associated production. 

As is well known, the CP-even can decay into W^W~, while the CP-odd cannot. The off-shell 
W^* from cascade decays must have the same charge as one of the two Ws from The leptonic 
decays of these two like-sign VFs result in the like-sign dilepton processes: 


pp^ i^i^Xrbbjj, (54) 

pp-> W^*H°A^ W^*W+W-bb i^e^Xrbbjj. (55) 


From Fig. I^we read off BR(Tf‘^ —)• W~^W~) ss 0.3, thus the theoretical cross section for this signal is about 
1.2 fb from and 0.6 fb from H^A^. The signals from resemble very much the background 

ttW^ since the like-sign FFs are usually well separated. The signals from H^A^ can be different from the 
background tiW^, but it is negligible compared to the same signals from fb) that we discussed 

in Sec. V A| If there are two Ff°s, they can both decay into W~^W~, inducing a like-sign trilepton signal. 


pp 


w^*w+w-w+w- 


(56) 


Although this signal is very distinct and background free, its production cross section is too tiny (0.037 fb) 
to be feasible. 


C. Signals for neutral scalars 

The Higgs pair production is an important process in SM as it can be used to probe the electroweak 
symmetry breaking sector. The dominant production mechanism at a hadron collider is through gluon 
fusion |[76l47^ . At LHC14, its cross section is about 18 fb at leading order, and is enhanced to about 
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33 fb by the next-to-leading-order QCD corrections 11781 . Recently, the next-to-next-to-leading-order QCD 
corrections are available, resulting in a further increase of the cross section to about 40 fb Il80l . We refer to 
Refs, llsn [821 for a more detailed discussion on the theoretical status of Higgs couplings at LHC. Physics 
beyond SM could alter Higgs pair production dramatically. The neutral Higgs pair production has been 
investigated in a variety of models including, for instance, MSSM IfTTl . NMSSM lIMlIMl . 2HDM ll85l - IMl . 
extended colored scalars Il85l [89l l90l . extended colored fermions Il85l 1911 l92l . Little Higgs |[93l . Higgs 
portal i94l |95l and composite Higgs models i9^ WT\ . There are also model independent approaches based 
on effective operators ||98]41001 . All these studies concentrate on the gluon fusion mechanism, since it is 
the dominant contribution. 

The associated production via the Drell-Yan process through an s-channel exchange of the Z 

boson has received less attention in previous studies due to several reasons. The single production of 
neutral scalars proceeds through gluon fusion and thus dominates in a large portion of parameter space. The 
associated production of neutral scalars through a Z boson exchange is much smaller than the gluon fusion 
channel in most models that we mentioned. And the charged scalars usually have more distinct features 
that may help to detect them at a collider. The circumstance could be modified significantly in the negative 
scenario of type II seesaw. Since the triplet vev is constrained by the p parameter to be small, the 
mixing between the doublet and triplet scalars and thus the couplings of /A^ to quarks are suppressed. 
As a consequence, the single and pair or associated production of the neutral heavy scalars through gluon 
fusion is also suppressed lIlOll . Since in the negative scenario more charged scalars are heavier than less 
charged ones, this makes associated production of neutral scalars comparable to those involving charged 
scalars. The matter becomes even more interesting when all neutral scalars are nearly degenerate, i.e.. 


K, Mj^o ~ Mfi- We therefore have chosen M^jo = 130 GeV in Eq. (30 1 as our benchmark point, 
and in addition we shall assume va = 10“^ GeV. The cross section for the associated H^AP production 
through the Z boson at LHC 14 is 350 fb at leading order, which is about 10 times as large as the SM Higgs 
pair production. 

Some recent studies 11811110241 1061 have demonstrated the potential of Higgs pair production in the signal 
channels 6677 , bhW~^W~, bbT~^T~ , and bbbb. With L = 3000 fb“^ at LHC14, the trilinear coupling of 
the Higgs bosons could be measured at an accuracy of ~ 40% II107II . With a much higher production cross 
section, we have carried out the simulation of H^A^ production at LHC14@300. To illustrate the potential, 
we shall consider the channel 0 


PP 


Z* 


0 .o 


H^A 


bbr'^i 


( 57 ) 


^ A systematic analysis on probing associated production of neutral scalars at LHC will be presented elsewhere [mu. 
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Signals like 6677 , hbW^W~ are also interesting, and will be studied in the future. An important part of 
this analysis depends on the ability to reconstruct the b pair and r pair. In our simulation, the 6 -tagging 
efficiency is assumed to be 0.7, and the misidentification rate of a c- and light-jet as a 6 -jet is taken to be 
0.1 and 0.01 respectively. We consider hadronic decays of r for its tagging, and assume an efficiency of 0.8 
with a negligible fake rate. 

The main SM backgrounds are: 


pp—)• bbZ/^*/h —)■ bbT~^T~ , 

(58) 

pp—> bbW~^W~ —^ bbT^h'rT~UT, 

(59) 

pp— ^ Zh — )• bbT~^T~ . 

(60) 


The irreducible background comes from bbT^T~, where the r-pair originates from the decay of Zj^* jh. 
Since the hadronic decays of r always have neutrinos, we also include the background bbW~^W~, which 
contributes to the bbr^VrT~ final sfafe. The bbW^W~ background dominantly stems from the ti pro¬ 
duction with subsequent decays t —)• bW and W —)• tVt-. In our simulation, we also include the associated 
Zh production with subsequent decays h ^ bb and Z —)• r+r~ or vice versa. 


Cuts 


bbr^T 

bbW+W- 

Zh 


S/B 

S/^/S + B 

Basic Cuts 

6797 

1529101 

952937 

2617 

4.31 

0.00274 

(N 

II 

cf 

II 

297 

20058 

34959 

102 

1.26 

0.00539 

ARbb < 2.0, < 2.0, > 1.5 

147 

2285 

3088 

35.4 

1.97 

0.0272 

> 80 GeV, > 40 GeV 

108 

537 

763 

21.3 

2.86 

0.0817 

Py' > 70 GeV, > 30 GeV 

68.5 

221 

93.7 

11.4 

3.45 

0.210 

Ht > 250 GeV 

44.0 

103 

30.3 

5.82 

3.25 

0.317 

95 GeV< Mrr < 135 GeV 

16.0 

7.94 

7.42 

0.426 

2.84 

1.02 

110GeV< Mbb < 150 GeV 

12.41 

2.65 

2.86 

0.284 

2.91 

2.14 


TABLE VII. Survival numbers of signal bbr'^r and its backgrounds, statistical significance S/\/S + B, and signal 
to background ratio S/B after imposing each cut sequentially at LHC14@300. 


We first employ some basic cuts for the selection of events: 


> 20 GeV, ARjj > 0.4, \7]j\ < 2.5. 


(61) 


Here the jets include both 6 - and r-jets. The numbers of the signal and background events after imposing 


the cuts are summarized in Table VII The last two columns in the table show the statistical significance 
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S/^/S + B and the signal to background ratio S/B. Then we apply the cut: 

Nb = 2, Nr = 2, (62) 

so that all four jets are tagged in this channel. At this stage, the background is about 200 times larger 
than the signal, since both and bbW~^W~ receive large QCD contributions although the fully 

electroweak Zh background is about 1/3 of the signal. To further suppress the background, we shall 
apply several kinematical cuts. In Fig. [^we show the distributions in the eight kinematical variables: 
the separations Ai2f,b,rT,fer. the total hadronic transverse momentum Ht, the leading (next-to-leading) 6-jet 
transverse momentum and the leading (next-to-leading) r-jet transverse momentum 

Because the two signal 6-jets come from decays of heavy neutral scalars, their separation ARhb is usually 
small with a clear peak around 1.4. In contrast, the two 6-jets in the bbT^T~ and bbW~^W~ background 
respectively are mainly from the gluon fusion and the ti decays, and thus tend to be separated. The same 
features are shared by the separation ARrr for the signal and the bbT~^T~ and bbW~^W~ background. For 
the Zh background, the separations ARbb and ARrr have a relatively large peak coverage from 1.5 to 3. 
Since the 6- and r-jets in the bbW^W~ background both come from t/t, their separation ARbr has more 
chances to be smaller than from other sources, as is seen in Fig. [T^ These differences between the signal 
and backgrounds suggest the following comprehensive cuts: 

ARbb < 2.0, ARrr < 2.0, ARbr > 1-5. (63) 

The distributions in the transverse momentum of various jets also provide useful information. First, 
both and in bbT~^T~ background are usually much softer, and the leading 6-jet transverse momentum 
p^ in the other two backgrounds is also small compared to the signal. Second, the r-jets in the signal 

are harder than in the backgrounds, since they come from the decays of heavier neutral scalars. We thus 
apply the following cuts on the transverse momentum of various jets: 

> 80 GeV, > 40 GeV, p^l > 70 GeV, pT^ > 30 GeV. (64) 

Since (anti-)neutrinos carry away part of energy in the hadronic decays of r, the r-jets in the final states 
will be softer than the 6-jets. Therefore, the pT cuts applied on the r-jets are 10 GeV smaller than on the 
6 -jets. Till this point, the QCD backgrounds are still about five times larger than the signal. Since the total 
hadronic transverse momentum Ht in the signal is usually larger than in the background, we further apply 
the cut: 


Ht > 250 GeV. 


(65) 
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FIG. 17. Distributions of various separations AR{,b^rT,bT and transverse momenta Ht, Px ^ p'r 

imposing basic cuts. See text for notations. 
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FIG. 18. Distributions of invariant masses in signal bbr^r at LHC14@300. Left panel; Mtt before cut ( 661 ; right 


panel: M^i, after cut ( 661 . 


After this cut, the survival bbW^W~ events become smaller than the signal, while the bbT~^T~ background 
is still more than twice as large as the signal. 

Now we consider to utilize the invariant mass distributions to further suppress the background. The left 


panel of Fig. 18 shows the invariant mass M^-r of the r pair. The bbT~^T~ background has a peak around 
M-rr = 70 GeV, which is the shifted peak of the Z boson pole at Mz = 90 GeV, and also has a sharp edge 
at Mz- In contrast, M^t in the signal has a much fatter peak. We therefore impose a cut on Mtt- 


95 GeV < Mrr < 135 GeV. 


( 66 ) 


This cut turns out to be very efficient. The reason that we apply a cut on Mrr first is due to the fact that Mrr 
is not only fatter than but also shifted from the original value of Mfjo/Ao according to our analysis. 


In the right panel of Fig. 18 we show the invariant mass of the bb pair. A sharp peak around 130 GeV is 
evident, so that the mass of H^/A^ can be reconstructed in this channel, although the CP nature of these 
scalars remains to be discerned. We then select the mass window of M^i,: 


no GeV < Mbb < 150 GeV. 


(67) 


After all these cuts, about 12 signal events survive atLHC14@300, with a statistical significance of 2.91 and 
a signal to background ratio of 2.14. So this channel is expected to be discovered with higher luminosity. 


VI. DISCUSSIONS AND CONCLUSIONS 

The new scalars in type-II seesaw are generically nondegenerate, when cascade decays between scalars 
proceed with the radiation of a VF boson. This can modify significantly their signatures at LHC even if the 
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mass splitting is not large so that the radiated W boson is far off-shell. In this work, we have presented 
a systematic and comprehensive analysis on this issue, by simulating all potentially interesting signals to 
the detector level and elaborately designing specific cuts to purify the signals. Our results may serve as a 
reference that encourages more detailed simulations by experimental collaborations. 

The positive scenario with doubly-charged scalars being the lightest is stringently constrained by LHC 
limits, and the scalars must be heavier than about 400 GeV. In contrast, the negative scenario is less 
restricted by current experiments. A much lighter spectrum is allowed in this case, and in particular, the 
new neutral scalars can be nearly degenerate with the standard model Higgs boson. This makes it very 
appealing from the point of view of collider phenomenology. 

Our main results are summarized as follows. In the positive scenario, the four-lepton signal is still 
the most promising channel to discover At LHC14@300, we can potentially probe up to a mass 

Mu±± ~ 600 GeV for NH and Mu±± ~ 700 GeV for IH. The five-lepton signal is crucial to deter¬ 
mine the spectrum of triplet scalars; for instance, a 5 cj significance can be reached wifh L = 76 fb“^ af 
Mfj±± = 400 GeV, Mjj± = 430 GeV for IH. The LNV signal from cascade decays of neufral scalars 
LfO/AO is heavily affecfed by fheir mass difference and fofal decay widfhs Tj^^o/^o due fo fhe 

significanf inferference befween fhem. For degenerafe and A^, our signal receives an enhancemenf fac¬ 
tor of fwo while fhe like-sign four-lepfon signal considered in fhe previous liferafure fends fo diminish. Buf 
ifs discovery af LHC 14 would require a large infegrafed luminosity due fo small cross section. 

For fhe negafive scenario, fhe associafed producfion of neufral scalars can give fhe same signals 

{bbT~^T~, bbW~^W~, and 6677 ) as fhe SM Higgs pair production, buf wifh a fen times as large cross 
section. We could reach 3a significance in ifs bbT^T~ signal channel af LHC14@300. The mosf promising 
signal is i^^Tbbbb from associafed producfion: a 5 cj significance can be reached af L = 109 fb“^ 

for LHC 14. Anofher signal from is i^^TbbT~^T~, which has a cleaner background and could be 

defecfed for L ~ 300 fb~^. The doubly-charged scalars being fhe heaviesf are relatively hard fo 
probe. We invesfigafed in defail fhe signal from associafed producfion, and found 

fhaf a 5a significance would require an infegrafed luminosify of L ~ 500 fb“^. 
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